Proposal of an experimental scheme for determination of penetration depth of 
transverse spin current by a nonlocal spin valve 
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We theoretically propose an experiment to determine the penetration depth of a transverse spin 
current using a nonlocal spin valve with three ferromagnetic (F) layers, where the Fi, F2, and F3 
layers act as the spin injector, detector, and absorber, respectively. We show that the penetration 
depth can be evaluated by measuring the dependence of the spin signal (magnetoresistance) on the 
thickness of the F3 layer. 

PACS numbers: 72.25.Ba, 72.25.Mk, 75.47.De 



I. INTRODUCTION 

There has been great deal of attention paid recently 
to the spin transport in nano-structured ferromagnetic 
materials because of their potential application to spin- 
tronics devices such as magnetic random access mem- 
ory and microwave oscillators. The giant and tunnel 
magnetoresistance effects^— and spin torque effect^ are 
key physical phenomena for the operation of these de- 
vices. The origin of these phenomena is the spin de- 
pendent electron transport, i.e., the spin current. The 
relaxation of the spin current in nonmagnetic (N) ma- 
terials has been investigated using, for example, a non- 
local spin valve system 8 -—. In ferromagnetic (F) mate- 
rials, on the other hand, we should distinguish the lon- 
gitudinal and the transverse spin currents, whose spin 
polarizations are parallel and perpendicular to the local 
magnetization, and which are the origins of the mag- 
netoresistance effect&ii and the spin torque effec t 6 ' d 8 , 
respectively. Compared to the extensive studies on the 
spin relaxation of the longitudinal spin curren t 3 ' 17 : 19 " ^!, 
there have been very few studies on the spin relaxation 
of the transverse spin current. 

The penetration depth of the transverse spin cur- 
rent (accumulation) is an important quantity character- 
izing its spatial spin relaxatior*22r~ . In our previous 
studies^r— , we proposed that the penetration depth was 
evaluated by spin pumping effect—; i.e., the creation of a 
pure spin current by ferromagnetic resonance (FMR), in 
the F/N/F trilayers. The point is that one F layer is in 
resonance and acts as the spin injector while the other F 
layer is out of resonance and acts as the spin absorber. 
Thus, by measuring the dependence of the FMR power 
spectrum of the injector on the thickness of the absorber, 
the penetration depth of the absorber can be evaluated. 
Because microfabrication is unnecessary, the FMR mea- 
surement is useful for investigating the spin relaxation. 
However, the choice of the F material is restricted to 
those that can to avoid the simultaneous resonance of the 
two ferromagnetic layers. Also, because the magnitude 
of the magnetization, as well as the FMR frequency, of 
the absorber increases with increasing its thickness, the 
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FIG. 1: Schematic view of the nonlocal spin valve with three 
ferromagnetic (F) layers. The black arrow in each F layer 
represents the magnetization. The red arrow indicates the 
spin of the conduction electron. The electric current flows 
from the Fi to the N layer while the pure spin current flows 
to the F2 and F3 layers. 



FMR spectrum of the injector and absorber sometimes 
overlap2I. These factors make it difficult to evaluate the 
penetration depth of the transverse spin current. 

In this paper, we propose an alternative method to 
determine the penetration depth based on the nonlocal 
geometry shown in Fig. [TJ The three ferromagnetic lay- 
ers are attached to the nonmagnetic layer, where the Fi, 
F2 and F3 layers act as the injector, detector, and ab- 
sorber of the spin current. By fabricating the three F 
layers in different shapes, we can control the directions 
of the magnetization in each layer independently; thus, 
a noncollinear alignment of the magnetizations can be 
achieved. For example, in Fig. [1] the F3 layer is assumed 
to be a cylinder with zero in-plane shape anisotropy; 
thus, the magnetization can rotate in the plane, while 
the magnetizations of the Fi and F2 layers are parallel 
to the stripe due to the shape anisotropy. Let us assume 
that the magnetization of the F3 layer is perpendicular 
to those of the Fi and F3 layers. Then, the amount 
of the magnetoresistance measured in the F2 layer de- 
pends on the absorption of the transverse spin current in 
the F3 layer. Thus, by measuring the dependence of the 
magnetoresistance on the thickness of the F3 layer, its 
penetration depth can be evaluated. There is no restric- 
tion in the choice of the materials and the change of the 
magnitude of the magnetization does not affect the mea- 
surement, which are advantages compared to the FMR 
method. 

The paper is organized as follows. In Sec. UH we discuss 
how to calculate the spin current at each F/N interface 
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FIG. 2: (a) The top and (b) side views of the system. 

by solving the diffusion equation of the spin accumula- 
tion. Then, the amount of the magnetoresistance in the 
conventional system is quantitatively estimated in Sec. 
IIII In this section, we also show that the penetration 
depth of the transverse spin current can be evaluated by 
nonlocal geometry. Section IIVI provides a summary of 
the work. The details of the calculations are shown in 
the Appendix. 



II. SPIN CURRENT AND SPIN 
ACCUMULATION 

In this section, we show the details of the calculation 
of the spin current and spin accumulation in the F and 
N layers, which are required to evaluate the magnetore- 
sistance effect. 

The details of the system we consider are schematically 
shown in Figs. [2] (a) and (b). The x and z axes are 
taken to be parallel and normal to the ./V layer, whose 
origins are located at the Fi/N interface. The thickness 
and width of the F& layer are denoted as d,F k and wp k 
while those of the N layer are c?n and wn . The distance 
between the Fi (F2) and the F3 layers are denoted as 
L\ (L2), respectively. The unit vector pointing in the 
direction of the magnetization in the F& layer {k — 1,2,3) 
is denoted as m^. The electric current, /, flows from 
the Fi to the N layer. By passing through the Fi layer, 
the conduction electrons induce spin accumulation inside 
and the interface of the Fi and N layers. These spin 
accumulations are diffused in the N layer, creating a pure 
spin current into the F2 and F3 layers. Then, the spin 
accumulations are created in these F layers. It should 
be noted that the following formula is applicable to an 
arbitrary alignment of the magnetizations, (mi, m2, 1x13). 

As shown in Sec. IIII1 the magnetoresistance effect is 



determined by the spin current at the F/N interface. Ac- 
cording to the spin dependent formula 1 -, the spin current 

at the F/N interface (into N, see lJ fe/N (k = 1,2,3) in 
Fig. [UJb)) is given by 



jF/N. 



y^ffF/Nm - (Mf - Mn) + ^7f /N-f F/N 



9% /N) m x (/*N x m) -gjj^ /N) /i N x 



t^ N jinx (juf xm)+t^/iF xm 



(1) 



where 7 F / N is the electric current from the F to N layer, 
9f /n = .9^ + g is the sum of the spin-up and spin-down 
conductance, 7 = (g^ — g^)/{g^ + g^) is the spin po- 
larization of the interface conductance 3 , is the real 

(imaginary) part of the mixing conductance^, and t^ju 
is the real (imaginary) part of the transmission mixing 
conductanc e 22 ' 24 ' 25 . The spin accumulation in the N and 
F layers are denoted as /^n and /^f, respectively. The 
first term (oc m) in Eq. (JTJ) describes the longitudinal 
spin transport (|| m), while the other terms describe the 
transverse spin transport (_!_ m). The terms proportional 
to g'^ describe the transverse spin injection from the N 
to F layer, while the terms proportional to describe 
the opposite flow of the spin current. In the zero pene- 
tration depth limit, the spin accumulation in the F layer 
is parallel to m, and the last two terms (oc t^) in Eq. ((T|) 
can be neglected, which is assumed in Ref. 18 . Although 
Refill assumes a spatially uniform spin accumulation, it 
has been shown that circuit theory is applicable to the 
diffusive syste m 28 ' 29 . 

The spin accumulation in the N layer obeys the diffu- 
sion equation 3 



dx 2 



\2 

At. 



(2) 



where An is the spin diffusion length of the N layer. The 
solution of is expressed as a linear combination of 
e ±x/\ N ^ rpj^ g pj n accumulation is related to the spin 
current by 



T 



dx 2e 2 



-A*n, 



(3) 



where 5n = wnc^n and <7n are the cross sectional area 
and the conductivity of the N layer, respectively. 

The longitudinal spin accumulation in the F layer, 
(m • /xp)ni, also obeys the diffusion equation, 



d 2 l_ 1 L 



(4) 



and its solution is expressed as a linear combination of 
e ±z / Ap , where Af is the spin diffusion length of the F 
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layer. The longitudinal spin accumulation is related to 
the longitudinal spin current by 



d HSf 
"dz 



i i 



(5) 



where S-p — iufWn is the cross sectional area of the F 
layer—. The electrochemical potential- and the conduc- 
tivity of the spin-cr electron are denoted as /i F and <7 F , 
respectively. The spin polarization of the conductivity 



is given by /3 = (cr p - <7 F )/(cr F 



<7 F ). The resistivity is 

defined as pp = 1/ctf = l/(cp + Op)- 

The transverse spin accumulation in the F layer, fip — 
m x (fi F x m) obeys the following diffusion equation 31 : 



dz 2 



T 
/1 F 



f 

A?' 



x 2 

F(T) 



T 
Mp. 



(6) 



(D F + Dp)h/(2J) relates 



The coherence length Xj = w ^ F -r j^ f j 
to the spin-dependent diffusion constant _D F and the ex- 
change interaction J between the conduction and local 
electrons^!. The spin polarization of the diffusion con- 



D, 



O/Pf 



For sim- 



stant is given by /3' = (D F 
plicity, we assume that /3' = /3. Then the transverse spin 
diffusion length A F ( T ) is given by A F (t) — Af/ yl — /3^21. 
The solution of Eq. ((6} is expressed as a linear combina- 
tion of e ±z / £ + and e _±z ^~, where 



f 



\ 2 \ 2 ' 

A F(T) a j 



(7) 



The penetration depth of the transverse spin current, At, 
is defined as 



-!- = Re 
At 



(8) 



The transverse spin current in the F layer relates to the 
spin accumulation as follows: 



m x 



x m = — 



d hSpap^ T 
"MF' 



dz 2e 2 



(9) 



where < = [4/(1 + /T) + a F /(l ~ p')]/2. 

We assume that the spin current is continuous at each 
F/N interface, and that the electric current is constant. 
Thus, the spin currents in the N layer near the interfaces, 
denoted as (i = 1 — 6) in Fig. [2fb), should satisfy 
the relation I^ 2 "" 1 + lf fc/N = I^a* (fc = 1, 2, 3). We also 
assume that at the ends of the N layer, the spin current is 
zero. Using these boundary conditions, we can solve the 
diffusion equations of fip and /xn, where the spin current 
at the F/N interface can be chosen as the integral con- 
stant of the diffusion equation of the spin accumulation. 
Then, the spin accumulations on the right hand side of 
Eq. ([1]) can be expressed in terms of the spin currents 
at each F/N interface. Thus, Is^ N can be obtained by 
solving their simultaneous equation; see Appendix A. 
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FIG. 3: The dependence of the magneto-resistance on the 
relative angle between mi and 1x13. The magnetization of the 
F2 layer is parallel to mi. 

III. MAGNETORESISTANCE 



In this section, we show the details of the calculation 
of the magnetoresistance. Also, we show that the pen- 
etration depth of the transverse spin current of the F3 
layer can be estimated by measuring the dependence of 
the magnetoresistance on the thickness of the F3 layer. 

By using the spin currents at each F/N interface ob- 
tained in the previous section, the magnetoresistance 
measured in the F 2 layer is calculated as R = R^ + 
R^ 2 \ where R^ and R^ represent the contributions 
to the resistance from the spin dependent transports in- 
side the F2 layer and the interface at the F2/N layers, 
respectively 3 -* 3 ^. The explicit forms of R^ and R^ are 
given by 



Pf 2 Sf 2 , 
2el 1 

2e 2 gF 2 



m 2 • /x F2 (z = 
d F2 



-d F2 ) - m 2 • /2 F2 (z = 0)] 



tanh 



2Aj 



(10) 



RW = - 



7f 2 /n 



2el 

7f 2 /n^Sf 



SF 2 m 2 • fi F2 (z = 0) - S N m 2 ■ fijsi{x = L)] 
1 2e 



2e 2 g F2 
2e 2 #N 



tanh(c?F 2 /AF 2 ) hi 
m 2 • I 



m 2 



. jF 2 /N 



(11) 



where gF/SV = h(l — /3 F )/(2e 2 pFAp) and ^n/Sn = 
h/{2e 2 p^X^), respectively. It should be noted that the 
above formula reproduces the results of Takahashi and 
Maekawa— by neglecting the F3 layer and assuming that 
mi = ±ni2. 

Figure [3] shows an example of the calculation result us- 
ing Eqs. (fTU|) and (fTTj) . The magnetizations of the Fi and 
F2 layers are assumed to be parallel (mi — 1112) while 
the magnetization of the F3 layer changes its direction 
from 1113 = mi to 1113 = —mi. The material parameters 
of the Fi, F2 and F3 layers are assumed to be identi- 
cal, for simplicity, and are taken to be p-p — 220 fmm, 
Af = 5.0 nm, Aj = 2.8 nm, f3p = 0.35, pn = 120 nm, 
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FIG. 4: The dependences of the magnetoresistance on the 
thickness of the F3 layer for the various penetration depth. 
The magnetizations of the F2 and F3 layers are assumed to 
be rri2 || mi and 013 _L mi, respectively. The dashed line is 
obtained by neglecting the transverse spin accumulation. 

A N = 1300 nm, R F/N = hS F /(e 2 g F/N ) = 275 A Onm 2 , 

7F/N 
,-2 



°- 70 > 9^ Fm /S F = 38.0 nm- 2 , gJ^ m /S F = 1.0 



nm t;.{ F/N) /S F = 4.u nm -, r i(F/N) /o F 
d F = 30 nm, c?n = 150 nm, w Fl — w F2 — 120 nm, 
10f 3 = 150 nm, wn = 150nm, L\ — 400 nm, and L2 = 
400 mn£.Z±22£±&M.. Then, the penetration depth is ob- 
tained as At = 3.5 nm. The resistance area, Eqs. (ITU)) and 
(fTTj) . is converted to resistance by using S F2 — w F2 Wn, 
as measured in the experimenter—. 

As shown in Fig. [3J the resistance, as well as the 
amounts of spin accumulation, of the parallel and an- 
tiparallel (1113 = ±mi) alignments are larger than those 
of the perpendicular alignment (ma _L mi), because of 
the fast relaxation of the transverse spin accumulation 
compared to that of the longitudinal one and because of 
the continuity of the spin accumulation. 

In Fig. rjj we show the dependences of the magnetore- 
sistance on the thickness of the F3 layer, d F3 for various 
penetration depths. The values of the parameters except 
d Fa and Aj are the same as those in Fig. [3J The pen- 
etration depth is changed by changing the value of Aj, 
where (Aj,A t ) =(0.7,1.0), (1.1,1.5), (1.5,2.0), (1.9,2.5), 
(2.3,3.0), (2.8,3.5), (3.5,4.0) nm, respectively. The mag- 
netization alignment is assumed to be 1112 = mi and 
ni3 _L mi. The magnetoresistance decreases as the thick- 
ness of the F3 layer increases because the amount of the 
spin accumulation created around the F2 decreases due 
to the spin absorption in the F3 layer. The magnetore- 
sistance is constant for a sufficiently large thickness of 
d Fs 3> At- The reduction of the magnetoresistance is 
fast for a small At because of the fast absorption of the 
spin accumulation. We also show the magnetoresistance 
obtained by neglecting the transverse spin accumulation. 
In this case, since the spin absorption in the F3 layer 
occurs only at the F3/N interface, and the amount of 
the spin absorption is independent of its thickness, the 
magnetoresistance is also independent of the thickness. 
In other words, the thickness dependence of the magne- 
toresistance reflects the relaxation of the transverse spin 
current in the F3 layer. Thus, the penetration depth of 
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n 
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4.0 nm 



the F3 layer can be evaluated. 



IV. SUMMARY 

In summary, we show a theoretical formula to calcu- 
late the magnetoresistance in a nonlocal spin valve with 
three ferromagnetic layers. We show that the penetration 
depth of the transverse spin current can be evaluated by 
measuring the dependence of the magnetoresistance on 
the thickness of the spin absorber. 
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Appendix A: Equations to Determine Spin Currents 

Here we show the details of the calculation of the spin 
currents at each F/N interface to evaluate Eqs. (|10p and 
PT|) . First, since the magnetoresistance depends on the 
relative directions of the magnetizations, we introduce 
three sets of unit vectors, (e^, ejty, m/.) (k — 1,2,3), 
which satisfy e^ x x &ky — m^. The rotational transfor- 
mation from (b2xj e 2y, 1*12) ^0 (bi x , eiy, mi) is character- 
ized by the rotation matrix Ri, which depends on the two 
angle parameters (#1, tpi) and is given by 



f cos 6*i cos (fi cos 9\ sin ipi — sin 1 ' 
Ri = ( — sinyji cost^i 
^ sin Q\ cos ipi sin B\ sin ipi cos 61 



(12) 



Similarly, we introduce the rotation matrix R3 which rep- 
resents the relative direction of {e 3x , e 3y , m 3 ) with re- 
spect to (e2x, e 2y, 1M2) and depends on (63, ips). Although 
the choice of the direction of the transverse unit vectors, 
ekx and e^y, is somewhat arbitrary, the final results is 
independent of these choice. The relation of the longitu- 
dinal and transverse spin currents between the different 
F layers can be given by (9k, <Pfe)^- 

The general solution of the diffusion equation of the 
longitudinal spin accumulation is given by 



4tt 



m • fi F 



g F sinh(dp/AF) 



7(1) 



r(2) 



2e 
H/3 F 
' ~2e 



I cosh 



cosh ( — 
Af 



(13) 



where Is} and li 2 J are the longitudinal spin current at 
z = and z — d F , respectively, flowing in the positive z 
direction. In the present study, as shown in Fig. [5] (b), 



I, 



(i) 



-m ■ lf^ N and I? 



0. The electric current / is 



positive for the electron flow in the positive z direction, 



5 



and is nonzero only in the Fi layer. The general solution 
of the spin accumulation in the N layer can be obtained 
by replacing the quantities of the F layer in Eq. ([T3|) with 
those of the N layer (/? in the N layer is zero). 

Similarly, the general solution of the transverse spin 
accumulation is given by 



e x ■ fi F = 47r 



(14) 



-ihlg + /i4 y 1} ~ ihlif + hl$] , (15) 



where fj(z) (J = 1 — 4) are given by 



h(z)=Re 



f 2 (z) = ilm 



cosh[(z — d F )/£^ 



g t smh(d F /£ + ) 
cosh[(z — d F )/£^ 



g t sinh(d F /£ + ) 



(16) 



(17) 



Appendix B: Explicit Form of Coefficient Matrix M 

Here we show the explicit form of the non-zero com- 
ponents of the coefficient matrix M , 

Mll = lH ^-^VrJffF./N (l-7g l/N )gFi/N 



2g Fl tanh(d Fl /A Fl ) 



2#n 



Mi 4 = - 



i 1 ~ T^ 1 /n)Sfi/n 



2.9n 



(23) 



I2,2 = ivi3,3 : 



"^(Fj/N) 1 ™ 



i'(Fi/N) 



Re 



3t(P0 tanh(d Fl /4-i) 



fft(Fi) tanh(d Fl /^ Fl ). 

(24) 



h(z) 



-Re 



cosh(z/f+) 



fi{z) = -ilm 



gt sinh(d F /£ + ) 
cosh(z/£ + ) 



g t sinh(d F /£ + ) 



(18) 



(19) 



Here fft /S F = [h/(4e 2 p F £ + )][(l+(3)/(l+f3') + (l-/3)/(l- 
iil and ijff are the transverse spin current at 

(2) (2) 

z = while I y s J and P a J are those at z — d F . In the 



present study, I. 



-e-r ■ I^ N and 7,„ 



while = = 0. 

By using the boundary conditions mentioned in the 
main text and the solutions of the spin accumulation 
above, we obtain the following simultaneous equations 
of the spin currents, 



r(2) 



-ej, ■ I 



F/N 



Ma = — s. 

2e 



(20) 



Here, M is the 18 th degree coefficient matrix, whose ex- 
plicit form is given in Appendix B. The 18 th de grcc vec- 
tor a consists of the spin currents at each interface, and 

, T Fi/N T Fi/N 

is given by ai = mi • Is 1 ' , a 2 = e lx ■ I s , a 3 = 
e ly ■ lf l/N , a 4 = mi • I^ 2 , a 5 = e lx • I^ 2 , a 6 = • I^ 2 , 

T F 2 /N T F2/N T F2/N 

a 7 = m 2 • I s , a 8 = e 2x ■ I s , a 9 = e 2y ■ h , 
aio = m 2 • 1^, a n = e 2x • I**, a 12 = e 2y ■ I**, 

r F 3 /N 

5 ' *s j a 15 — c 3w ' ■ 

rN 6 _ „_ tN 6 



T F 3 /N tFs/N r r 

ai3 = m 3 ■ l s , a i4 = e 3x ■ l s , a i5 = e 3y ■ l s 
ai 6 = m 3 -I^ 6 , a i7 = e 3a ,-I^ 6 , a i8 = e 3j ,-I^ 6 . The source 
term of the spin current and spin accumulation is given 
by the 18 th degree vector s, whose non-zero component 
is only Si, 



, i 1 ~ 7 Fi/n )5f 1 /n ( d Pl 
si = 7Fi/N + P Fl ; tanh 



2 ffF! 



2A 



Fi 



(21) 



Then, by numerically calculating the inverse of M, the 
spin currents at each F/N interface, a, can be obtained. 



'2,3 : 



_ gj(Fx/N) ti Tm 
-M 3 ,2 -— * r(F i/N) Im 



.+t4- Dp 

^(Fi/N)^ 



1 



g t ( Fl ) tanh(d Fl /^ Fl )_ 
1 



fft(Fi) tanh(rf Fl /£ Fl ) 



M 2 ,5 = M 3 , 6 = - 



M 2 , 6 = -M 



3,5 = 



g r(Fi/N) 
5N 



g i(F!/N) 
5N 



M 4 ,i = M 5 ,2 = M 6i3 = 1, 



M44 = M 55 = M 



6,6 : 



1 



tanh(Li/AN) 



M 4 ,13 = -M 4 ,16 : 



sinh(Li/AN) 



1 



M 4 ,15 = -M 4 ,18 = 



sinh(Li/Aj. 



1 

sinh(ii/AN 



M 5i i 4 = -M 5! i7: 



1 



(25) 
(26) 

(27) 
(28) 
(29) 



(RxrV 1 ),,, (30) 



M 4;14 = -M 4 ,i7 = - 3-^ — ( R i R 3 X )3i. (3 1 ) 



■(RxR,- 1 )^, (32) 



M 5 ,i 3 = -M 5 , 16 = - . 1 . (RiRg 1 )^, (33) 
smh(Li/AN) 



sinhCLi/A,)^ 1 ^" 1 )-' 
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M 5 ,i5 = -M 5 ,i8 = -- 



1 



sinh(ii/AN) 



(R1R3 1 )i2i (35) 



Mm. 



18 : 



1 



sinh(L 2 /A N ) 



(R 3 _1 )32, (49) 



M 6 ,i3 = -M 6i i 6 = 



sinh(ii/AN) 



(RiR3 1 )23, (36) 



M.,u = -M.,ir = - ail|h(£ 1 i/Aw) (RiR,- 1 ) a i, (37) 



M 6 ,i5 = -M 6i i 8 = - 



sinh(Li/AN) 



(RiR 3 " 1 ) 2 2, (38) 



111,16 = — 



sinh(L 2 /AN) 



(Rs 1 )^, 



Mn, 17 = 



sinh(L2/Af 



;( R 3 )ii) 



Mn, is— . —77 — x (R3 
sinh(L 2 /AN) 



(50) 
(51) 
(52) 



M7J = 1 + 



( 1 ~ ^ 2/n ).9f 2 /n 
2#f 2 tanh(rf F2 /A F2 ) : 



(39) 



M ^' = M(W ^ a ' (53) 



M 7 , 



10 : 



(! ~ 7p 2 /N)ffF 2 /N 

2.9N 



(40) 



Mi2,i7= . * ■ r (R 3 " 1 )2i, (54) 
sinn(iv2/ANj 



M 8>8 = M 9j9 =l+i^ F2/N) Re 



+ *i(F 2 /N)I m 



3t(F 2 )tanh(d F2 /^F 2 ) 
1 



5t(F 2 ) tanh(cfr 2 /4v,) 



(41) 



Ml2,18= . 77 J R 3 ' 1 )22, (55) 

smh(L 2 /A N ) 



_ (1 - 7g 3/N )gF 3 /N . 
2g N smh(Li/A N ) 



M 8i9 = -M 9;8 = -i^ F2/N) Im 



+ *i(F 2 /N)^ C 



1 



5t(F 2 )tanh(d F2 /^F 2 
1 

5t(F 2 )tanh(c? F2 /^ F2 ) 



M 8) n = M 9) i2 = 



M 8 ,l 2 = -M 9 ,n: 



#r(F 2 /N) 
.9N 



g i(F 2 /N) 



(42) 



(43) 



(44) 



M l : i -.= ^"j^ ) ^ (R3Rr 1 )3i, (57) 



M 



13,6 = 



2g N sinh(Li/A N ) 



^ ~ 7F 3 /N)gF 3 /N 

2g N sinh(ii/A N ) 



(R 3 Rr 1 )32, (58) 



Ml3,13 = l + 



( X ~ 7 F3 / n ).9f 3 /n (1 - 7 F3 /n)3f 3 /n 



+ 



2<7 F 3 tanh(d F3 /A F .,) 2gN tanh(Li/AN) ' 

(59) 



M 



13,16 : 



( 1 ~ 7 F3 / N )gF 3 /N 

2.g N tanh(Li/A N ) ' 



(60) 



Ml0,7 = Mii,8 = Mi2,9 = 



tanh(L 2 /AN) ' 



Mi ,io = Mii,ii = Mi2,i2 



1 + 



tanh(L 2 /AN) 



Mm,i6— . —77 7T n (R3 X )33, 
sinh(L 2 /AN) 



(45) 



(46) 



(47) 



M 



M 



14.4 — 



14,5 : 



M 



10,17 = 



sinh(L 2 /A N ) 



(R^si. (48) 



* 3/N) -(R 3 Rr 1 )i3+ ffijf* , (R 3 Rr 1 )23, 



g N sinh(Li/A N ) 



g N sinh(Li/A N ) 



(61) 



5rt(ip3/N) -(R 3 Rr 1 )n+ , g &f> , (R 3 Rr 1 ) 2 i 



g N sinh(Li/A N ) 



#n sinh(Li/A N ) 



(62) 



Ml4,6 = 



5rt(ip3/N) -(R 3 Rr 1 )i2+ g fe; r /N > ^ (R3Rr 1 ) 22 , 



3 N sinh(L 1 /A N ) 



g N sinh(L 1 /A N ) 



(63) 



7 



fjl ff r(F 3 /N) , n 

Ml5,4— ■ u/r /\ \ l K 3 K l J23- 



#i(F 3 /N) 



3 N sinh( J L 1 /A N ) ' g N sinh(L 1 /A N ) 



M 15 , 6 = J RaRr^i- 



?i(F 3 /N) 



5 N sinh(Li/A N ) ' ,g N sinh(Li/A N ) 



T4- 

g r (F 3 /N) rp R _n 

IVI 15,6 — _ — — 1_ / t- — ~ — 7l K 3 K l )22- 



#i(F 3 /N) 



5 N sinh(L 1 /A N ) g N sinh(Li/A N ) 



(R3R1 )l3i 
(64) 

(RsRj^ 1 )!!, 

(65) 

(R3Rr 1 )i2i 

(66) 



Mi4,i4 = Mi5,i 5 =l + 



#r(F 3 /N) 



g N tanh(Li/A N ) 



^ r(F 3 /N) 



+^k/N) im 



1 



5t(F 3 )tanh(d F3 /^F 3 ) 
1 



5t(F 3 )tanh(c? F 3/^F 3 ) 



(67) 



M 



1 



17.6 : 



sinh(Li/AN) 



(RsRr 1 )^, (76) 



Ml8,4= ■ , ,1 , (R 3 Rr 1 )23, (77) 

sinh(Li/ANj 



Ml8,5 = " 



1 



sinh(Li/AN) 



Ml8,6 = 



sinh(Li/AN) 



(R 3 Rr 1 ) 2 i, (78) 



(R 3 Rr 1 ) 22 , (79) 



Mi6,7 = -Mi 6 ,io = -- 



1 



sinh(Li/AN) 



(Rs)33, (80) 



Ml6,8 = -Mi 6! n = - 



sinh(Li/Aj> 



-(R 3 ) 3 i, (81) 



Mi4,15 = -Mi 5 .i4 : 



y i(F 3 /N) 



<?n tanh(Li/AN) 



r r(F 3 /N) im 



+ *i(F 3 /N)^- C 



5f t(F3) tanh(dF 3 /^F 3 ) 
1 



5t(F 3 )tanh(d F3 /^F 3 ) 



Ml4,17 = Mi5,i 8 = 



Mi4,18 = -Mi 5 ,17 = 



1 



1 U 

'r(F 3 /N) 



SfNtanh(ii/A N ) ' 



#i(F 3 /N) 



5 N tanh(Li/A N )' 



M 



M 



M 



M 



16,5 : 



16,6 : 



17,4 : 



sinh(Li/A 



1 

sinh(Li/AN) 



sinh(Ll/A^ 



"(R3R1 1 )32, 



sinh(Li/AN) 



17,5 = 



sinh(Li/AN) 



(R3R1 1 )i3i 



(R3R1 )iij 



(68) 



(69) 



(70) 



Ml6,4 = -^- 7r - T -T( R 3 R r 1 )33, (71) 



(RaRr'jsi, (72) 



(73) 
(74) 
(75) 



Mi6,9 = -Mi 6 12 = 



Mi 7 7 = — M17. 



10 = 



Mi8,7 = -Mi 8 10 = 



Ml8,8 = -Mi8,ll = 



Mi8,9 = -Mi 8 12 = 



sinh(Li/AN) 



(R 3 )32, (82) 



sinh(Ll/A^ 



;(R 3 )i3, (83) 



Mi7,8 = -Mi 7! n = - 



Mi7,9 = -Mi 7! i2 = - 



sinh(Li/AN) 



(Rs)n, (84) 



sinh(Li/Aj> 



-(R 3 )i2, (85) 



1 

sinh(Li/AN) 



(Rs)23, (86) 



sinh(Li/Ai> 



-(R 3 ) 2 i, (87) 



sinh(Li/AN) 



Mi6.13 = Mi7 i i4 = Mi 8 15 = 



M 16,16 = Mi7 i i7 = Mi 8> i 8 : 



tanh(Li/AN) ' 



(R 3 ) 22 , (88) 



(89) 



1 1 
■ + - 



tanh(Li/AN) tanh(L 2 /AN) 
(90) " 
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